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While it is clear that in some kinematic regime QCD can be described by an effective 
(as opposed to fundamental) string theory, it is not at all clear how this string theory 
should be. The 'natural' candidate, the bosonic string, leads to amplitudes with the 
usual problems related to the existence of the tachyon, absence of the adequate Adler 
zero, and massless vector particles, not to mention the conformal anomaly. The super- 
symmetric version does not really solve most of these problems. For a long time it has 
been believed that the solution of at least some of these difficulties is associated to a 
proper identification of the vacuum, but this program has remained elusive. We show in 
this work how the first three problems can be avoided, by using a sigma model approach 
where excitations above the correct (chirally non-invariant) QCD vacuum are identified. 
At the leading order in a derivative expansion we recover the non-linear sigma model 
of pion interactions. At the next-to-leading order the 0(p 4 ) Lagrangian of Gasser and 
Leutwyler is obtained, with values for the coefficients that match the observed values. 
We also discuss some issues related to the conformal anomaly. 
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1. Introduction 

It seems that a paper with this title should necessarily start with a list of the 
reasons that make_us believe that there should be an effective description of QCD 
in terms of stringsEJ. 

The more commonly cited arguments are the dominance of planar diagrams in 
the large N limita 'filling in' a surface (interpreted as the world-sheet of a string), 
the expansion in terms of surfaces built out of plaquettes in strong-coupling lattice 
QCDcl, and the success of Regge phenomenology!!!, which can ultimately be under- 
stood in terms of string theory ideas (although, as we will discuss a little bit later, 
the actual Regge theory that corresponds to QCD cannot be derived, at present, 
from any known string theory). 

To these we could add two more reasons. One_is the appearance in string theory 
of the universal (at long distances) Liischer termcl. The static inter-quark potential 
provided by the string V(r) = ar + c gets modified by quantum fluctuations by a 
Coulomb-like piece — 7r/12r. While for the heavy quark bound state calculations it 
is perturbative QCD (complemented with some non-perturbative corrections) that 
is applicable, and not the string regime, the Liischer term is quite useful. Finally, 
and in a completely different context, namely that of deep inelastic scattering, the 
evolution of the parton distribution down to low values of Q 2 (around (2 GeV) 2 ) 
leadsB to a low x behavior for the structure functions of the form x~ 1A7 , while 
Regge theory predicts x~ l , in striking good agreement. 

In the last years the possibility of describing some non-abelian gauge theories in 
terms of string theories has become a reality thanks to the AdS/CFT relationshipO. 
Unfortunately, the gauge theories that one can describe, or rather solve, in this way 
possess an unrealistically large number of supersymmetries. Breaking this large 
amount of supersymmetry to the realistic case of N = while decoupling the un- 
wanted degrees of freedom from long-distance physics has proven to be a formidable 
task and one in which not much success has been attained yet. 

It should be clear to the reader from the very beginning that we are not ad- 
dressing in this paper these deep and fundamental issues. Yet, it is clear that the 
mentioned theoretical developments give credence to the idea that there must be 
some kinematic regime where the simplest bosonic string, obtained after integrating 
out all the remaining degrees of freedom made heavy by the breaking of supersym- 
metry, must provide an approximately valid description. 

In fact it is highly questionable where one should expect a string description 
to be exactly valid for real, non-supersymmetric, four-dimensional, asymptotically 
free QCD, as a string does not seem to be the natural language to understand high- 
energy processes in deep- inelastic scattering where the point-like structure of quarks 
and gluons is apparent. We should probably be less ambitious and satisfy ourselves 
with an effective description. We subscribe this point of view and think of strings 
as effective theories and not worry at all about their mathematical consistency as 
fundamental objects. 
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It is surprising that even with this limited scope all known string theories fail to 
provide a description of low energy QCD. To see how this comes about let us recall 
the original Veneziano amplitudeO. The original motivation for this ad hoc formula 
was to provide a description of hadrons that manifestly exhibited duality. While thii 
evidence for duality was rather weaktl, at the time this amplitude sparked Nambuli2l 
to suggest that it should originate in bosonic string theory. 

After decorating the Veneziano amplitude with the appropriate Chan-Patonliil 
factors^, it is supposed to describe the scattering amplitude of four pions 

A{ir a ir b -> ttV*) ~ Tr(T a T b T c T d )A(s , t) + non cyclic permutations, (1) 

1 (— a{s) — a(t)) 

where a(s) = 1 + a's is the Regge trajectory. The Regge trajectory slope is deter- 
mined, after a fit, to be a' ~ 0.9 GeV~ 2 El We immediately recognize that there are 
poles in the s-channel whenever a's = n — 1. Thus a tachyon is present for n = 
and we see the first of the problems alluded to in the abstract appearing. Further- 
more the second (J = 1) state (which we should identify with the rho particle) is 
massless. This is the second one of the problems we alluded to. 

Shortly afterward, supersymmetry was introduced in string theorya, but unfor- 
tunately most of the difficulties with the bosonic string are still present in super- 
symmetric string theory when one applies it to the present context. The relevant 
amplitude for our purposes is the four tachyon amplitude |^| 

A( ~_ T(l - a p (s))T(l - a p (t)) 
T(l - a p (s) - a p (t)) 

Now ct p (s) = 1 + a's is the so-called p-trajectory which in Regge theory describes 
exchanges of particles with positive G-parity. This is the Lovelace-Shapiro ampli- 
tude and it contains no tachyonic poles. Could it then be a candidate to describe 
pion scattering? The answer is no. It does not have the appropriate Adler zero, 
i.e. the property that at s = t = the pion scattering amplitude vanishes. And 
although it has no tachyons in the intermediate channel, unfortunately, the vertex 
which is supposed to represent the pion is a tachyon itself. Furthermore, the first 
resonance in the p-trajectory (supposed to represent the physical p particle ) has 
zero masiL All in all, one is at the same dead end. 

A fidlj to this problem consists in arbitrarily changing the intercept in a p {s). 
Indeed, if we write a p {s) = | + a's and use Eq. (^), the resulting amplitude has 
poles in the s or t-channel when a's = n + h. It has no tachyons and the first 



a Due to difficulties with unitarity only orthogonal groups or certain representations of unitary 
groups can be introduced in this way. |— ■ 

In supersymmetric theories one usually performs the GSO projection! 1 !, projecting out the 
tachyon. However one may choose not to do so and compute the four-tachyon amplitude, supposed 
to describe scattering of spinless particles in the bosonic (Neveu-Schwarz) sector. 
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pole, which is massive, can be identified with the p particle. Furthermore, the 
previous amplitude has the right Adler zero. Expanding this amplitude in powers of 
s and t and comparing with the amplitudes obtained from an effective Lagrangian, 
Polyakov and VereshaginEl found that L\ = ^L 2 , L 2 = g^p-ln2, L 3 — 
— 2L2, which numerically turn out to be quite acceptable valuesQ Unfortunately, 
no known string theory leads to such a fix in the intercept^. 

It is believed that the ultimate reasonfpr the presence of a tachyon in the spec- 
trum lies in a wrong choice of the vacuun£3. Since the choice of the spin zero vertex 
operator, V{k) —: expikx :, is based on the Lorentz properties alone, it is the same 
both for scalar and pseudoscalars and, accordingly, both scalars and pseudoscalars 
have tachyonic poles in the s-channcl on account of parity conservation. The situa- 
tion is thus parallel to the one in multicomponent <f> A theory with V((f>) = — /i 2 +A</> 4 , 
where perturbing around = gives negative m 2 values for all components. It is 
natural to assume that the amplitudes obtained through the use of the canonical 
vertex operators correspond to (unphysical) amplitudes for excitations perturbed 
around the wrong, unphysical vacuum. These ideas have in fact been around for a 
long time, but no one appears to have implemented them in any practical way. 

In a previous paperc3 by some of the present authors a possible line of attack 
was suggested. The idea is the following. Given that with the present string theory 
technology it is impossible to find the 'true' vacuum, let us assume that the vacuum 
is non-perturbative in string creation and annihilation operators and that it actually 
corresponds to the true QCD vacuum. In this case, the relevant (lightest) degrees of 
freedom are the ones emerging after the spontaneous breaking of chiral symmetry. 
In the physical vacuum of QCD there is a clear distinction between scalars (sigma 
particle) and pseudoscalars (pions). The sigma particle is just another hadronic 
resonance with a mass in the 1 GeV region and without any specially important role 
to play. The massless pseudoscalars, on the other hand, can be collected in a unitary 
matrix U(x) which under chiral transformations belonging to SU(3)l x SU(3)r 
transforms as U(x) — > U'(x) = LU{x)R^ and describes excitations around the non- 
perturbative vacuum. From the string point of view U{x) is nothing but a bunch 
of couplings involving the string variable x. Our goal is to find a consistent string 
propagation in this non-perturbative background. 

A crucial property of string theory is, undoubtedly, conformal invariance. This 
amounts to demanding that the theory is independent of the specific conformal 
factor chosen to describe the two-dimensional world sheet. While this is a desirable 



c The relation L\ =y-U2L2 = — I/4L3 was established earlier in bosonization models and in the 
chiral quark model£3 by means of a derivative expansion of quark determinant. However at that 
time its possible connection with a string description of QCD was not recognized. 
d It appears possible to modify to some extent the Regge trajectory for low values of J, M 2 |-apd 
hence the intercept by adding additional terms in the string action, rigidity, for instance ! 1 ^H 20 ! 
or imposing constraints on the string motions£J. However spectroscopy suggest that while such 
terms might exist they are probably of limited influence. Moreover the conformal invariance of a 
modified string becomes questionable. 
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and, as a matter of fact, a crucial property of fundamental strings, it might not be 
necessarily so for the QCD string (if we look at the QCD string with a magnifying 
glass we shall eventually see quarks and gluons, not the string itself!) . However in 
the limit of large N c the QCD hadronic amplitudes are saturated by gluon fishnet 
diagrams among which any higher order samples give, in principle, a comparable 
contribution. Thereby it appears plausible that there should be no dependence 
on a specific choice of configuration describing a dominating gluon propagation. 
Therefore the hadronic string action should obey reparameterization invariance of 
diagram surfaces and conformal invariance. Since conformal invariance must hold 
when perturbing the string around any vacuum, perturbative or not, we have a 
powerful tool at our disposal, namely to couple the chiral field U (x) to the string 
degrees of freedom and demand conformal invariance. Exactly in the same way as 
Einstein equations are derived from string theorycj by requiring the vanishing of 
the beta function for an external metric G fjLl ,{x). 

This approach was followed in E3 with mixed success. The lowest order chiral 
Lagrangian, or rather the equations of motion coming from it, were obtained and 
the whole approach seemed consistent. Unfortunately, the 0(p A ) coefficients were 
identically zero. Experimentally these are numbers of order 10 -3 . We have now 
found out the reason for this failure. It turns out that the action considered in §1 
is only a particular case of the general case, and one that does not allow a proper 
treatment of the unitarity condition of the external source U{x). When the general 
case is considered, unitarity and locality can be implemented consistently — at least 
at the order we have worked — and, as a consequence of the new terms in the action, 
new divergences and hence new contributions to the beta function for U{x) appear. 
These lead to 0(p 4 ) terms in the effective Lagrangian. These are the results we 
would like to report here. 



2. Coupling pions to the QCD string 

The hadronic string in the conformal gauge is described by the conformal field 
theory action in the four dimensional Euclidean space-time ^] 

W s tr = T^7 / d 2+t v (-) diXpdiXn, (4) 



Ana 1 J \fi 
where for e = one takes 

x n — x fi ( T , v) ; — °o < T < oo, < a < oo; i — r,a |Lt = l,...,4. 

The conformal factor ip(r, a) is introduced to restore the conformal invariance in 
2 + e dimensions, and this is the only way it enters the theory. The Regge trajectory 

c The direct evidences for dimensions of the string world shcct P^lH (d=2) and of the target spaceFl 
(D=4) were found from the analysis of meson state densities at high energies as compared to the 
Hagedorn-like experimental growth of the latter ones. 
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slope (related to the inverse string tension) is known to be universal a 1 ~ 0.9 GeV -2 

EE 

We would like to couple in a chiral invariant manner the matrix in flavor space 
U (x) containing the meson fields to the string degrees of freedom while preserving 
general covariance in the two dimensional coordinates and conformal invariancc 
under local scale transformations of the two-dimensional metric tensor. The equa- 
tions of motion for the U(x) field will then be obtained from the condition that the 
quantum theory must be conformally invariant, i.e. the (3 functional for the U(x) 
couplings must vanish. 

Since the string variable x does not contain any flavor dependence, we have 
to invent a way to couple it to the background U(x) variable. We introduce two 
dimensionless Grassmann variables ('quarks'), or rather several families of them, 
living on the boundary of the string sheet: ip l{t) , i> r{t) which transform in the 
fundamental representation of the light flavor group (SU(2) in the present paper). 
A local hermitean action Sb = J drLf is introduced on the boundary <j = to 
describe the interaction with background chiral fields U(x(t)) = exp(i7r(cc)// 7r ) 
where the normalization scale /„• ~ 93MeV, the weak pion decay constant, is 
introduced to relate the field n(x) to a 7r-meson one. 

The boundary Lagrangian is chosen to be reparameterization invariant and in 
its simplest minimal form reads 



L 




f = h, (i> L u{i - z)ip R - ip L u{i + z)Tp R 

+i> R U+(l + z*)^ L - i> R U + {\ - z*)tl> L ) , (5) 

herein and further on a dot implies a r derivative: ip = dtp /dr. In order to arrive to 
a number of field redefinitions have been made. It is impossible to simplify 
any further. Details are given in Appendix A. 
A further restriction is obtained by requiring CP invariance. There are two 
CP-like transformations. The first one is 

U U+, IpL «-> ipR- (6) 

The above Lagrangian is CP symmetric for z — —z* — ia. The second one is 

U^U+, ^ L ^U+^p L , ip R ^Uip R . (7) 

Under this transformation the Lagrangian becomes invariant only for z — 0. We 
interpret the first CP transformation as the physical one and the one which one 
should require of a Lagrangian describing strong interactions. 

The above coupling may appear surprising at first and somewhat ad hoc. To 
see that this is not so, let us expand the non-linear field U(x), i.e. U(x) ~ 1 + 
in(x)/f n + ... and retain the first two terms. The first term just gives rise to a 
^-function propagator which eventually leads to the familiar ordering in the usual 
string amplitudes t\ < <2 < The second term just provides (after integrating the 
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fermions out) the usual (tachyonic!) vertex. In short, if we ignore the non-linearities 
in the theory we are back to the usual difficulties. 

It is easy to see that the previous action is invariant under general coordinate 
transformations of the two dimensional world sheet. The fermion action is auto- 
matically conformally invariant, because it does not contain the two dimensional 
world sheet metric tensor since it can be written as a line integral. 

InE3 the cases z = ±1 were considered. None of them is a valid one on symmetry 
grounds. 

3. Diagrammar 

Now we expand U(x(t)) around a constant background xq and look for the poten- 
tially divergent one particle irreducible diagrams (1PI). We classify them according 
to the number of loops. Each additional loop comes with a power of a! . 

We expand the function U(x) in powers of the string coordinate field x^(t) = 
XOfi + %ft (t) around a constant xq which is the translational zero mode of the string 

U(x) = U(x ) + x |1 (r)0 f4 ?7(xo) + -x^Xu^d^duUixo) + ... 

= U(x )+V(x). (8) 

One can find a resemblance to the familiar derivative expansion of chiral perturba- 
tion theory. Indeed perturbation theory in the operators (H) makes sense as a low 
momentum expansion which is presumably valid up to momenta approaching to 
the mass of the first massive resonance {p meson etc.). In the present case a' is the 
dimensional parameter normalizing the above expansion. 
The free fermion propagator is 

(^r(t)Mt')) = U- 1 (x )9(t - t'). (9) 

If we impose CP symmetry then 

(Mt)Mt')) = 0>r(t)iPl(t')) 1 = U(x o )0(t - t'), (10) 

for unitary chiral fields U(x). 

The free boson propagator projected on the boundary is 

(^(rK(r')) -^A(r-r'), A(t -f r') = A(0) ~ -, d T A(r -> r') - 0, (11) 

e 

the latter results hold in dimensional rcgularization (see below). 

In order to make contact between dimensional regularization, a short-distance 
cut-off (which we shall later use) and Regge phenomenology we need to unam- 
biguously fix the normalization of the string propagator. This can be inferred 
from th e definition of the kernel of the N-point tachyon amplitude for the open 
stringEl The Veneziano amplitude corresponds to the insertion of vertex operators 
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: exp(ifc^ x^{rj)) : on the boundary of the string. After resolution of the Gaussian 
integral one obtains for the kernel of the generalized beta-function 



exp( 



V Si* 



3>l 



2a'fe 0) fe (i) 



which unambiguously prescribes 

A fa _ T i) = -2a'ln(|rj - r/|/i). 



(12) 



(13) 



The \i dependence does not show up in (|12| ) due to energy-momentum conservation. 

Keeping in mind this definition let us determine the string propagator in dimen- 
sional regularization, restricted on the boundary. First we calculate the momentum 
integral in 2 + e dimensions 



A £ (r) 



4W ( - 



d 2+e k exp(ik r) 



) J (2tt) 2 + £ k 2 



2a' 



e V ip 



O(e) 



(14) 



A dimensionally regularized propagator properly normalized to reproduce (|13|) can 
be constructed by subtracting from ( llif ) its value at t/x = 1 where ( |l3| ) should 
vanish 



Ae(r) 



c— 



)( 




— e 




1 












i|r/x| 


+ 0(e) 









Therefrom one finds unambiguously the relation 



A(0) 



-2a' 



1 



+ C + In <p 



(15) 



■0(e) = A £ -2a'ln<£>, (16) 



where in the spirit of dimensional regularization we have assumed that e < and 
hence the first term in (|l5|) vanishes at r = 0. 

The two-fermion, iV-boson vertex operators are generated by the expansion 
(||) and they appear with an extra sign (— ) = i 2 following the definition for the 
generating functional Zf, = (exp(iSb)) and Eq. (|^). In particular, for the L — > R 
transition one has 



V=~\ ((1 - z)V(x)d T + (1 + z)d T [V(x) ...]). 



(17) 
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and for the R L transition the hermitean conjugated vertex V + appears. The 
corresponding Feynman rule for the 2-fermion, N-boson vertex of the L — > R tran- 
sition is 

V N = -^d A 6(A -B)d^... d^ N U(x ) [(1 - z)S(A - n) • • • 6{A - r N )+ 

+ (1 + z)S(t 1 -B)---6(t n -B)}(18) 

where A, B are proper-time values for the left- and right-handed fermions and 
n, . . . , rjv are proper-time values for the boson field-string variables. 

From this point on it is quite straightforward to proceed with the renormal- 
ization process. We shall determine the counterterms required to make the beta 
functional for the coupling U(x) vanish up to the two loop level. In spite of the rel- 
ative complexity of the Feynman rules, the fact that we are working with a boundary 
field theory is crucial in making the calculation manageable. In fact most diagrams 
can be determined by simply playing with integration by parts and using basic 
properties of the Dirac delta function. Yet the renormalization is quite non-trivial 
and the ultraviolet structure of the counterterms is surprisingly quite complex. It 
is thanks to this complexity that non-zero values for the 0(p ) coefficients can be 
obtained. In fact we believe that some of the results presented in this work can have 
some bearing on more general discussions involving fundamental strings too. 

In what follows we retain not only the singular parts of one-loop diagrams but 
also the finite ones as they will be necessary to construct two-loop diagrams. 



4. Renormalization of the fermion propagator at one loop 

To avoid clutter the main body of the paper we have relegated the detailed deriva- 
tion of the different Feynman diagrams to the appendixes. Since the present calcu- 
lation is somewhat non-standard, we provide the technical details there. 

Using the set of Feynman rules described in the previous section one arrives to 
the following result for the divergent part of the propagator (Appendix B), 

3 + z 2 



OiA-B^AW- 1 



9 -d^uu- 1 d^u\u-\ 



(19) 



This divergence is eliminated by introducing an appropriate counterterm U 
U + 5U 



SU = A(0) 



2 p 



j-dJJU-^U 



= 0. 



(20) 



Conformal symmetry is restored (the beta- function is zero) if the above contribution 
vanishes. 

Let us find out for which value of z this variation of U is compatible with its 
unitarity. 



5(UU+) = U ■ SU + + SU ■ U + = 0. 



(21) 
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A simple calculation shows that this takes place for z = ±i. For other values of 
z eq.(pp|) entails d^UU^ILU = which has only a trivial constant solution in 
Euclidean space-time. In c3 this value of z was not considered and thus unitarity 
of U was not properly taken into account. 

When z = ±i and before the unitarity constraints are imposed the local classical 
action which has ( pp| ) as equation of motion is 



For other values of z ^ ±l;±z the related local action is unknown. The above 
Lagrangian is of course the well known non-linear sigma model which is commonly 
employed to describe pion interactions. 

We have thus succeeded in finding the action induced by the QCD string. It has 
all the required properties of locality, chiral symmetry and proper low momentum 
behavior (Adler zero). Furthermore, it describes massless pions. f^, the overall 
normalization scale, cannot be predicted from these arguments. 

To this point we have been quite successful in our program, but of course no real 
predictability has been achieved yet. Indeed we knew the form of this action from 
general principles, even though it is nice to see that things work out consistently. 
We have to turn to the 0(p 4 ) effective Lagrangian to get non-universal results. 

5. Renormalization of the vertices at one loop 

In order to proceed to a two loop calculation we shall need in addition to the 
counterterms for the one-loop propagator (which we just got) the counterterms for 
the vertex with two fermion lines and one and two boson lines, a; M , respectively. We 
also have to check whether the minimal Lagrangian (^J) is sufficient to renormalize 
also the vertices. It turns out that, in fact, it is not. 

Let us obtain the divergences for vertices with external boson lines. We introduce 
an external background boson field x^ to describe vertices with several boson legs 
and split x M = x^ + r)^. The free propagator for the fluctuation field rj^ coincides 
with the one for x^. 

Then the total divergence in the vertex with two fermions and one boson line is 
(see Appendix C) 



W&) = Ijl / d 4 xtr [d^UdJJ- 1 + d^U+d^U+r 1 } . 




(22) 



9(A - BJ-A^JI/- 1 + z) [-d^{d 2 U) + 2d v UU-%d„U 

+ (1 + z)d^d v UU- 1 d v U - -(1 + z)(3 - z)d I/ UU- 1 d lll UU- 1 d v U 
+x ll {B){\ - z) [-d^{d 2 U) + (1 - z)d v UU- 1 d ll d v U + 2d„d v UU- 1 d v U 
-i(l - z)(3 + z)d v UU- 1 d li UU- 1 d v U jt/- 1 



(23) 
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Let us now amputate the fermion legs and evaluate the pure vertex divergences. 
The amputation rules are 

Xfi (A)9(A -B) = - J dTd T 6(A - t)x^t)6(t - B), 

Xp{B)6{A -B)= J dr9(A - t^t^^t - B). (24) 



Then the divergent part (|23|) can be reproduced by the following operator in the 
Lagrangian 



(^l<5> {1 Hr-4>l<S> (2) ^r) +h.c, $^ 2 ) = x^t)^ 2 \ (25) 



i 

2 

where the vertex matrices ' can be rearranged as follows 



<&« = A(0) |(1 - z)d, (^-d 2 u - ^^-d v uu- l d u u 
+zd u UU~ 1 d ti UU- 1 d u U) } 

= (i - z)d ft (su) - 0„ 

^ = {l + z)d^(5U) + ^. (26) 

The terms proportional to derivatives of 8U are automatically eliminated by the 
redefinition of U that one performs to renormalize the one- loop propagator (and it of 
course vanishes if the equations of motion are imposed). But the part proportional 
to <j)^ remains and to absorb these divergences new counterterms are required. 
Evidently, these terms come out of the following terms in the Lagrangian: 



AL dm . = -A(0)(1 - z 2 )ip L \ -—d v UU-^d v U —d v UXJ- Y d v XJ 

+zd l ,UU- 1 UU~ 1 d u U^ ip R + h.c. (27) 

Therefore the counterterms required to eliminate the additional divergences for the 
vertex with one boson and two fermion lines can be parameterized with three bare 
constants g\ , g 2 and (73, which are real if the CP symmetry for z = —z* holds 

AL bare = i(l - z 2 )i, L - zg 2 )d u UU- 1 d v U - ( 9l + zg 2 )d u UU- 1 d u U 

+2zg 3 d u UU- 1 UU~ 1 d u U^ i> R + h.c. (28) 
Renormalization is accomplished by redefining the couplings gi in the following way 

9l = g i>r - A(0). (29) 

The constants Qi j- are finite, but in principle scheme dependent, and from ( |l6| ) it 
follows that a logarithmic dependence of the bare couplings on the conformal factor 
ip is introduced along the renormalization process. The counterterms are of higher 
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dimensionality than the original Lagrangian (^) and therefore the couplings gi are 
of dimension M~ 2 . Since (p|) was actually the most general coupling permitted 
by the symmetries of the model, in particular conformal invariance, one is lead to 
the conclusion that conformal symmetry is broken, already at tree level, by these 
couplings, unless they happen to vanish. Since they are dimensional, it is natural 
to normalize them by the only dimensional parameter, namely a' . 

Even if the new couplings are dimensional, it turns out that at the order we 
are computing, the trace of the energy-momentum tensor is still vanishing once 
the requirements of unitarity of U are taken into account (see Appendix D) and 
therefore conformal invariance is not broken at the order we are working. At higher 
orders in the a' expansion further counterterms may be however required in order to 
ensure conformal invariance perturbatively. We postpone a more detailed discussion 
to the final sections. 

One can introduce the running "effective" couplings 

9l + A e = g lir + 2a' In cp = gf, (30) 

One-loop conformal invariants are gf — g 2 and gf — g%- At any rate, the dependence 
of the new couplings on the Liouville mode is determined. 

In any case, the appearance of new vertices changes the fermion propagator due 
to the diagrams discussed in Appendix E. One obtains from such terms (which are 
of higher order in derivatives) the following contribution to the propagator 

9(A - B)i;A(0)(l - z 2 )U- x {2(g hr - z 2 g 2 , r )d p UU- 1 d p d p UIJ- 1 d p U 

-(1 + z)[gx, r + zg 2/r )d p UU- 1 d p UU- 1 d p d p U 
-(1 - z)(gi, r - zg 2tr )d p d^UU- 1 d p UU- 1 d p U 
+4z 2 g 3 , r d p UU- 1 d p UU- 1 d p UU- 1 d p U} U' 1 
= -9(A - B)A(0)U- 1 5 W UU- 1 (31) 

we shall denote this contribution by 8(A — B)d g . One should add the divergence 
contained in d g to the one-loop result, thereby modifying the U field renormalization 
and equations of motion 



SU = A(0) 



2 



= 0. (32) 



This is, in fact, the source of the much sought after 0(p 4 ) terms. 

We note thai this new contribution (^lj) is proportional to 1 — z 2 and it was 
thus absent in E3. At this point we have to ask whether such equation of motion 
may be derived from a local effective Lagrangian containing both dimension 2 and 
dimension 4 operators. This would then constitute the effective Lagrangian derived 
from the string model. However at this point this question is too premature to 
formulate. Two-loop diagrams generated from (|^) could certainly produce similar 
contributions and, as a matter of fact, so could new counterterms from diagrams 
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with two fermions and two boson lines, should they require an additional countert- 
erm. In fact it can be seen that the above equation of motion cannot be derived from 
a local Lagrangian involving the unitary matrix U . The requirements of locality and 
unitarity would force g^ r — 0, so this should not be the full answer. 

Before concluding this section and moving to the other contributions we have 
just mentioned, we calculate the renormalization of the vertex with two fermion 
and two boson lines, as this is also required as a counterterm. Let us summarize 
the divergent structure for these diagrams (see Appendix F) 

9(A - B)-^- 1 {x„(A)x v (A) [d^(-SU)(l + z) - 
+x^B)x v {B) [dpd v {-5U){l -z) + <t>»u] 

+A(0) fracl - z 2 2 / dr [-x fl (T)2 u (T)d p UU- 1 d fi UU- 1 d v d p U(l - z) 

J B 

+Xp(T)x u (T)d p dpUU- 1 d u UU- 1 d p U(l + z) 

Hxp(t)x u (t) - i ll {r)x v (r))d p d li UU- l a„d p U] } ET 1 , (33) 

where 



^ = A(0)- 



1 — z 1 4- z 

dpdpd ' " ' [ : ' ' " ' :) I :i : 



-d p dpd v uu- x d p u — —d p uu- l dMpU 



-zdpdpUU^dvdpU + zdpUU^d^dJJU^dpU 

+2z{dpUU- 1 dpUU- 1 d u dpU + dpdpUU^d^UU^dpU) 

~2zdpUU- 1 dpUU- 1 d u UU~ 1 dpU] (34) 

One can check that the terms encoded in (f)^ , together with the last contribution in 
Eq. (|3^) combine precisely as a second variation of the additional interaction vertices 
(1271) m coordinate fields. Therefore their renormalization is completely performed 



with the help of counterterms (28) and no additional counterterms or operators 
appear. 

It can also be seen that, in fact, any diagram with an arbitrary number of 
external boson lines and two fermion lines, i.e. any vertex of those generated by the 
perturbative expansion of (|J) is rendered finite by the previous counterterms. This 
completes the renormalization program at one loop. 



6. The fermion propagator at two loops 

The two-loop contributions to the fermion propagator can be obtained from one- 
loop diagrams with two external boson legs by joining the latter with a boson 
propagator. A factor 1/2 has to be added. One must include not only one-particle 
irreducible diagrams but also one-particle reducible ones and consider both diver- 
gent and finite parts. 

There are 10 two-loop diagrams which are listed in Appendix G. The divergences 
in the propagator at two-loops can be presented separated into five pieces 

9(A-B)[di + dn + dni + div + d v ]. (35) 
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The first and second piece contain the double divergence A 2 (0), the third, fourth 
and fifth pieces reveal only a single divergence A(0). Finite parts are irrelevant for 
the present discussion. 



Tabic 1. Chiral field structures proportional to the double divergence A 2 (0) 
appearing in the two loop contribution to the fermion propagator. 



CF structure dj djj Total 





P 2 P 2 




1 

8, 





1 

8_ 




P 2 P — p 




3+z 2 





3+z 2 




p — P 2 P 




8 2 

3+z 2 
8 , 





8 2 

3+z 2 
8 , 




p — p 2 — p 
p — p — PP 
PP — p — p 
PP — PP 




3+z 2 





3+z 2 






3+z 2 8 (3+z 2 ) 2 


(l-z 2 )(l + z) 2 


8 2 3 
ll + z + 5z 2 -z i 






8 32„ „ 
3+z 2 (3+z 2 ) 2 


32 

(l-z 2 )(l-z) 2 


16 

ll-z + 5z 2 +z 3 






8 , 32 
3+z 2 


32 



lf 2 

3+z 2 

9 8 2 


p 


— p — p - - 
p — PP — p 


■p 


3+z 2 . (3+z 2 ) 2 
8 "r , 32 
(3+z 2 f 

(3+z 1 ) 2 




(1-z 2 ) 2 


(3+z 2 )(7+z 2 ) 

32 2 
1 + z 2 


p 


— p — P-- 


p 


(l-z^)z 2 


(9-z 2 )(l + z 2 ) 


16 


8 


16 



The component di represents "the second variation" 
the one-loop divergence 



or one-loop divergence in 



3 + z 2 



SU = A(0) 



dl(5U)- 

[d^sup^d^u - d^uu- 1 su u-%u 

3 + z 2 



2 M 



-dJJU- x d*U 



■ dpUU-%(SUj\ \u-\ 
(36) 



Therefore it is renormalized away by the redefinition of the U field and vanishes 
when the equations of motion are imposed. The counterterms renormalizing U field 
yield the same expression but twice more andjif the opposite sign. Thus the result 
is — dj in correspondence with the results of E3 (for z = 1 they coincide). This is 
all that was obtained in that work. In particular no single-pole, A(0) appeared and 
therefore no new equations were obtained at the two loop level. Accordingly, the 
coefficients of the 0(p 4 ) coefficients were deemed to vanish. This will not be the 
case here. 

The second part represents the remaining terms of order A 2 (0) in two loop 
diagrams after subtraction of di and it reads 

1 



<I n = §A 2 (0)(1 - z 2 )V- x {2(1 - z^dpUU^d^dpUU^d^U 

-(1 + zfdpUU^d^UU^dpd^U - (1 - zfdpdfJJU^dpUU^d^U 
+Az 2 d p UlJ- 1 dp,UU- 1 d p UU- 1 dp,U} U' 1 . (37) 
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This term is identical, but of opposite sign, to the contributions generated by the 
one-loop counterterm in the vertex with two fermions and one boson line, after its 
insertion in a one- loop diagram (see Appendix E). 

To summarize, we show in Table 1 the distribution of A 2 (0) divergences between 
di and dn- Short-hand notations are used for the corresponding chiral field (CF) op- 
erators, for instance, p pp p corresponds to U~ 1 d^UU~ 1 d^dpUU~ 1 d p UU~ 1 . 

In dm we include those single-pole divergences, proportional to A(0), which 
are removed once the one-loop renormalization of U in the finite nonlocal part of 
fermion propagator at one loop (see, Eq. ( |B.2| )) is taken into account, that is when 
we replace U by U + S^U in the one- loop propagator 

i(l-2 a )A(0)A(AB)^ _1 [d»(5U)U-%U 

-d^UU' 1 SU U-^dyU + df.UU^d^SU)] U~ l . (38) 

Likewise in div we include the divergences that are eliminated when the additional 
counterterms in the one-boson vertices (those proportional to gi) are included in the 
finite part of the one- loop fermion propagator (the terms proportional to A(A,B) 



in Eq. flE.ll)). 

One can check that all terms in the two loop fermion propagator linear in A(0) 
and in A(A, B) belong either to dm or to d\y . We present them in Table 2. Thus 
one-loop renormalization removes dm and djy completely. 



Table 2. Summary of single divergences A(0) which are eliminated 
after the introduction of the one-loop counterterms. 

CF structure 



p fi* p 

P P PP 

PP P P 

PP PP 

P P P P 

P PP P 

P P P P 





dni 


dfv — 


Total 




l-z 2 





l-z 2 




1- 8 Z 2 
l 8 -z 2 






8 2 
l-z 2 

8 2 
l-z 2 


(1 


-z 2 )(3+z 2 ) 


(l_ z 2 )(1 + z) 2 


9 8 


(1 


-z 2 )(3+ 2 2 ) 


,16 „ 
(1- Z 2 )(1- Z ) 2 


9 8 

(l-z 2 )(l + z) 


(1- 


16 


z 2 )(3+z 2 ) 


16 




(l-z 2 ) 2 


8 



(l-z 2 )(3+z 2 ) 


(1 


- Z 2 )(3+. 2 ) 


16 

l-z 2 


(1- 


9 8 9 

z 2 )(3+z 2 ) 


8 9 9 

(l-z 2 )z 2 


9 2 9 

(l-z 2 )(3-z 2 ) 




8 


4 


8 



Some single-pole divergences remain however. Indeed, there are some diver- 
gences linear in A(0) which come from the double integral in the diagrams of 
Appendix G, HJ) and flG.10|) , 

J(A, B) = / d n [ 1 dr 2 d Tl A(n - r 2 )d T2 A(n - r 2 ) 

J B J B 
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A ~ B 2 

(It(A-B-t) A(t) . 



(39) 



In Appendix H this integral is calculated using two different regularizations. The 
divergence is found to be 



d v = cyA(O) [XJ- 1 d p UU- 1 d IJ ,UU- l dJJU- x dpUU- 



-U- x d n UU- 



l dpUU~ 



(40) 



z 2 ) 2 /8 — a'/2 for z = ±i. This term survives after adding all 



with cy — a'(\ 

the counterterms and together with (31) are the only new genuine divergences that 
can contribute to the beta function (single poles). It must therefore be added to 
the equation of motion at the next order in the a 1 expansion and modifies the term 
8^U , S^'U — > S^U + dv in a crucial manner; namely it opens the way to non 
zero solutions for the coupling constants gi and therefore for nonzero values for the 
Gasser-Leutwyler 0(p 4 ) coefficients. 



7. Local integrability and unitarity 

The equation of motion of 0(p 2 ), Eq. fl2C|), can be obtained from a local action of 
the Weinberg type fl22]), involving a unitary matrix U(x), only for z = ±i. If the 
corresponding terms with four derivatives that we have just found are to be derived 
from dimension-four operators in a local effective Lagrangian then certain relations 
are for sure to be required from the so far arbitrary constants gi >r . 

Such a Lagrangian has only two terms compatible with the chiral symmetry if 
we use the dimension- two equations of motion (|2^) , 

£< 4 ) = ~/> (K^UdpU^d^UdpU- 1 + ^Ud^dpUdpU- 1 + h.c.) . (41) 
The terms 

dlUdpU^dpUU- 1 , dlUd 2 p U-\ (dl) 2 UU-\ dpdpUd^dpU- 1 

which are in principle possible are reduced to the set (^lj) with the help of integration 
by parts in the action and of the dimension-two equations of motion (|2p| ) . 

Variation of the previous Lagrangian gives the following addition to the equa- 
tions of motion 

^ = -flu- 1 {2K X [dpd p uu- l dpUU- l dpU + d^uu^dpUu^dpdpU 

-df+UU^dpUU^df+UU^dpU - 2dpUU- 1 d p UU- 1 d p UU- 1 dpU 

+d, l uu- 1 d 2 p uu- 1 dpU] 

+K 2 [d^dpUU^dpUU^dpU + dpUU- 1 dpUU- 1 dpd p U+ 

2dpUU- 1 d f ,dpUU- 1 dpU 

+d 2 pUU- 1 dpUU- 1 dpU + d^UU^dpU^dp 



1, 2008 14:1 WSPC/Guidelines-IJMPA hadstr 



Hadronic string, conformal invariance and chiral symmetry 17 

-dJJU- 1 d p UU- 1 d p UU- l d ll U - Zd^UU^dpUU^d^UU^dpU 
-ZdJJU^dpUU^dpUU^dpU] } U~ l . (42) 



Table 3. Comparison between the coefficients of the different chiral field structures 
in the equations of motion derived from a local Lagrangian and from the condition 
of vanishing beta function, namely vanishing of the single pole divergences at the 
two loop level 



CF structure x - l a S r - d,g dy 





ft - 




P 


tip 




-2{2K X + K 2 ) 




z 2 )(l 


+ 2 )(9l,r + Z92,r) 







ftp 




-ft - 


-p 




-2(2K l +K 2 ) 




, 2 )(i 


- z )(9l,r - Z92,r) 







ft - 




Hp - 


-p 




-4K 2 






~9l, r + Z 2 g2,r) 





P 




P 


P 




P 


2[(1 - z 2 )K! + K 2 ] 







-cv 


P 




'P 


P 




P 


-2z 2 K 2 












n 




P 


ft 




P 


4[Ki + K 2 ] 




-1(1- 


- 2 2 )z 2 ff 3 ,r 


cv 



Let us now apply the 0(p 2 ) equations of motion to remove the Laplacian on 

chiral fields d 2 U . Then one obtains the set of coefficients for the various chiral field 

re- 
structures given in Table 3. These coefficients to be determined are then compared 

with the results obtained from the coefficients of the one-loop and two-loop single- 
pole divergences (see (|3l| ) and (f!o"|)). For z 2 = — 1 only one solution is possible, 
implying 

1 a 1 

K 2 = 0, K 1 =--C V = - — \ gi,r = -g2,r = -93,r = 4c V , (43) 

Thus, comparing Eq.f[4lJ) with the usual parameterization of the Gasser and 
Leutwyler Lagrangiar£3, 

L, = \l 2 = -\l, = -\k^ = ^. (44) 

For oi = 0.9 GeV_ 2 and /„ ~ 93 MeV it yields L 2 ~ 0.9 • 10~ 3 which is quite a 
satisfactory result! 3 .!!. 

So far we have paid no attention to the unitarity of U at the two-loop level. 
Does the variation implied by d g and dy respect the perturbative unitarity of U? 
Or does this lead to a new constraint eventually incompatible with the numerical 
previous numerical values? It turns out (see Appendix I) that if one accepts ar- 
bitrary real coefficients in the set of dimension-four vertices included in (^|) then 
the only solution compatible with the unitarity is given by the parameterization 
with constants K\ and K 2 . Thus the requirement to preserve unitarity under field 
renormalization is entirely equivalent to the local integrability condition, similarly 
to the case of dimension-two operators. This is a remarkable result that hints to 
the consistency of the whole procedure. 
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8. Conclusions: conformal invariance and all that 

Conformal invariance is a subtle issue. As we have discussed in detail in the intro- 
duction this is a necessary ingredient for the consistency of string theory and it is, 
according to general principles of field theory equivalent to requiring the vanishing 
of the beta functions. 

There are several sources of breaking of conformal invariance in the present 
model. We have been concerned with the equations of motion for the U(x) string 
field functional. The beta functional has been computed up to two loops, properly 
including the unitarity constraints. The vanishing of the single pole divergences at 
one loop leads to the familiar non-linear sigma model. At two loops the 0(p 4 ) terms 
of Gasser and Leutwyler appear. The numerical values turns out to be identically 
determined by the arguments of locality and unitarity and they therefore determine 
the 'renormalized' value of gi <r in terms of the constant cy. The latter being the co- 
efficient of a single loop divergence in an integral appears to be manifestly universal. 
The former could be in principle ambiguous since a change in the renormalization 
scheme (for instance by choosing the subtraction to a value different from fj/r = 1 
in the x-propagator) may shift a finite piece between A(0) and the renormalized 
value of the coupling. 

However, this ambiguity is of no relevance. The requirements of locality and 
unitarity seem to restrict the 'measurable' value of gi^ r to specific, well determined 
values. If one changes the renormalization scheme by a finite amount 5, the 'mea- 
surable' couplings will all change by the same amount to gi r = gi_ r + 6, but it is 
this last quantity the one that will then appear in the equations of motion and will 
be related to cy, which is universal, so this change of scheme is immaterial. 

Furthermore, the renormalized couplings <7, lT . are conformal factor independent, 
so the relations determined in the previous section hold in any conformal frame 
and so are the values for the 0(p 4 ) coefficients of the chiral Lagrangian. This is 
of course in full agreement with general considerations based in QCD and chiral 
Lagrangians (in the large N limit, that is, where the string picture is supposed to 

hold) -o 

In E3 a general, but tentative, argument (based in conformal invariance) was 
given that, if valid, would imply the vanishing of the 0(p ) coefficients to all orders 
in the inverse string tension. Here the argument simply fails because the additional 
couplings which are required to renormalize the model are dimensional. These cou- 
plings were not considered in E3 because the need for additional counterterms was 
not manifest along the renormalization process due to the (incorrect) choice of the 
value of z. 

Of course this immediately raises a new issue. The new counterterms seem to 
entail a breakdown of 'classical' conformal invariance. Indeed, introducing the cou- 
pling <?;, which are of dimension M~ 2 makes the action non-conformally invariant 
already at the classical level, in a way that is similar to the introduction of a cou- 
pling to the tachyon in the familiar bosonic string theory. Although the issue is 
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somewhat independent of the equations of motion for the matrix field U[x) we 
have to give an answer to this difficulty if conformal invariance is to hold. 

To demand conformal invariance in the string-mode sector, we have to guaran- 
tee the vanishing of the trace of the energy momentum tensor. This issue can be 
discussed in several ways, but let us proceed to determine the trace of the energy 
momentum tensor. On general grounds we can write 

6 = ]TA0 4 , (45) 

i 

where is the beta function for the i-th coupling and Oi the accompanying oper- 
ator. Then, when computing the trace of the energy momentum tensor, the part of 
the tensor bilinear in fcrmion fields leads to a contribution that contains operators 
of the form ( p8| ) with their couplings gi replaced by the appropriate beta func- 
tions. These beta functions contain a classical, purely engineering, part as well as a 
quantum one, which has been computed in the text, and which is the same for all 
couplings gi. In the string-mode sector their contribution is given by the averaging 
over the fermion vacuum and it vanishes according to the arguments presented in 
the Appendix D. There is no a fermion induced string action when the CP sym- 
metry (|^) holds. Thus the presence of conformally non-invariant interaction on the 
boundary does not affect conformal symmetry in the bulk. 

Throughout this paper we have systematically used conventional perturbation 
theory and the fermion loop expansion which however may be more cumbersome 
when one proceeds to the vertices with larger number of emitted bosons and to 
higher orders in the loop expansion. As an alternative, direct functional methods can 
be developed, an example of which is displayed in the Appendix J. This approach 
can be also exploited to show how to go ahead without fermions and so bypassing 
the potential problems related to breakdown of conformal anomaly since it allows 
a formulation where the 'quarks' are integrated out. Doing so apparently removes 
any conformal anomaly present in the boundary as we just discussed. 

On the other hand, at some point we do expect a genuine breakdown of confor- 
mal invariance. This is so because of the conformal anomaly that the bosonic string 
theory must necessarily exhibit. In fact, this breakdown (tantamount to introduc- 
ing a scale dependence) is quite welcome as we do not expect the simple string 
picture presented here to be even approximately true at very short distances (even 
within the large N limit). The naive bosonic string action used inthe present paper 
does not prevent large Euclidean world sheets from crumpling E3. It does not also 
describe correctly the high-temperature behavior of large N QCD El 

According to by now common loreEfl, the dependence on the conformal factor 
has to be understood in terms of the renormalization group flow. This dependence 
has not been worked out in detail, but it would be feasible to take it into account 
perturbatively in a way not too different from standard sigma model calculations 
in string theory by including dilaton, spin two, etc. backgrounds. The dilaton beta 
function precisely contains the d — 26 term characteristic of the anomaly. At some 
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order in the a' and topological expansion there would appear a coupling between the 
dilaton and the chiral fields through the exchange of x string variables. It would be 
interesting how the hadronic degrees of freedom can modify the anomaly condition, 
if at all. _____ 

A QCD induced string action may also include nonlocal ______ or, at least, 

higher-derivative __ vertices breaking manifestly conformal symmetry which help 
to make the strings smooth and supply the correct high-temperature asymptotics. 
However, as we are concerned here with the low-energy string properties we do 
not expect that the strategy and technique to derive the chiral field action needs 
any significant changes to be adjusted to a modified QCD string action. The only 
changes may come out from the short distance behavior of the modified string 
propagator. 

Finally we enumerate some of the simplifications and missing points of the 
approach undertaken in the present paper. 

1) From the beginning we have restricted ourselves to the SU(2) global flavor 
group. The reason is that only parity-even terms in the equations of motion can be 
revealed from the simple fermion Lagrangian (||) . As it is supposed that all relevant 
meson degrees of freedom are reproduced by the hadronic string one cannot expect 
the appearance of multi-fermion interaction which may effectively arise only due to 
heavy-mass reduction of glueballs and hybrids suppressed in the large N limit. Then 
the only way to obtain the parity-odd Wess-Zumino-Witten terms is to supplement 
one-dimensional fermions with spinor degrees of freedom, i.e., for instance, to add 
reparameterization invariant vertex 

AL f = Zftfx ll 1>, (46) 

with Dirac matrices 7 M and new dimensional constant £. This extension will be 
investigated elsewhere. 

2) One can easily add external electromagnetic fields and thereby calculate also the 
chiral constants Lq, Liq. However it is not yet clear how to introduce current quark 
masses consistently with the open string picture with ultra-relativistic quarks at 
their ends. Hence the other chiral constants need more efforts to be understood 
within the string approach if the boundary fermions were indeed associated with 
quarks. ______ 

3) Going backlo the possibility of higher derivative __ or nonlocal l____ or five- 
dimensional ____ action for the hadronic string, we find it interesting to explore the 
traces of deviations from the conformal string theory in low-energy chiral constants 
and phenomenology. Including the running of the 0(p 4 ) coefficients would however 
require the introduction of l/N corrections. We do not know if the string picture can 
be consistently implemented when subleading corrections are taken into account. 

4) As we have mentioned in this section — and emphasize here once more — the 
conformal anomaly can hopefully be taken into account consistently within the 
present approach, but he have not really attempted to do so in this present work. 
This is a major challenge that we leave for the future. 
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5) The value for the coefficients of 0(p 4 ) that we have found are, up to the usual 
factors of ir, etc. rational numbers. In other words, they are not the first terms 
in the Taylor expansion of a T function or similar transcendental functions. What 
would be the appropriate amplitude replacing Vcncziano formula then when the 
proper symmetries of the QCD vacuum are considered such as in the present case? 
We do not know. 

All things considered, we believe that the objectives we set in the introduction 
have been achieved. We have obtained the leading effective action for low energy 
QCD from two very simple requirements: chiral invariance and conformal symme- 
try and two simple ideas: perturbing about the true vacuum of QCD and using 
the simplest possible effective world-sheet action. The approach seems to be over- 
simplified, but has proven to be considerably robust. All cross checks and physical 
requirements are met. The outcome is very sensible from a phenomenological point 
of view. Perhaps more importantly, several interesting avenues are open for future 
exploration and a fully consistent approach may be developed. We also hope that 
the considerations presented in this work may be of use to the ample string com- 
munity since some of the techniques we have employed appear to be relevant in a 
wider context. 
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Appendix A. 

To construct the most general reparameterization- and chirally invariant action on 
the boundary of the string one can use the following set of operators bilinear in 
fermion variables and of the minimal dimension one 

tPlUtJjr, Tp L Uip R , tPrU+iPl, yL R U + ip L , 

i>LlpL, ^RlpR, 4> L UU + ll) L , ^rU+U^r. (A.l) 

Other vertices like tpLUtpR can be decomposed in a linear combination of basic 



vertices (A.l) after integration by parts in the action S& = JdrLf. The multi- 
fermion local interaction is suppressed in the leading large N approximation as we 
assume that all meson degrees of freedom relevant in this limit are reproduced by 
the hadronic string. Then a multi-fcrmion interaction may effectively arise only due 
to heavy- mass reduction of glueballs, hybrids and multiquark mesons suppressed in 
the large N limit. 

As the case of CP symmetric action is of the most importance to provide the 
conformal symmetry we restrict ourselves with the analysis of this action only. Thus 
the general CP invariant hermitean Lagrangian takes the form: 

Lf = bi> L Uvb R + b*i>LUip R + b*J R U + i> L + bvb R U + ^ L 

+ic [i^L^L + + id (^ l UU + vj l - ip R U+U^ , (A.2) 

where the constant b is complex whereas the constants c, d are real. As compared 
to the minimal Lagrangian (||) the general Lf contains three more real parameters, 
Im(6), c, d. Now let us show that by the local rotation of fermion variables preserving 
their chiral structure, 

tpL = aiipc + a^UtpR, vb L — a{rpL + ou^rJJ^ , 

^r = Mr + P2U + ^ l , $r=R$ r + I%$ l U ) (A.3) 
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Fig. 1. One- loop diagrams for the propagator. 



one can eliminate the redundant vertices and reduce the Lagrangian (|5|) to the 
minimal form. In order to prove it we transform the minimal Lagrangian (^|) with 



the help of rotation ( A.S ). The initial set of constants is b = (a + i)/2, c = d = 0. 
The final set of vertices will be CP invariant if the following conditions are fulfilled: 

ai/3f = alft, (A.4) 
a 2 p* 2 = «2/3 2 , (A.5) 
ai p* = a^Pi. (A.6) 

The first two constraints relate some phases, 

ai = |ai| exp(z0i), Pi = ±|/3i| exp(z0i), 

a 2 = |a 2 | exp(z0 2 ), #2 = ±\fo\ exp(# 2 ), (A. 7) 

whereas the third one eliminates one of the moduli of \cxj | and \(3j\. Three remaining 
moduli and the relative phase <j>\ — </> 2 turn out to be sufficient to fit three real 
constants Im(6), c, d, 



\ai\\fh\ = |a 3 ||/3i| = yd + ^±fl! = I C , 

Q 

COS(>i -<j) 2 ) = -, 

±| ai ||A| ± \*2\\H = \C ± jf}) = ^ Im W- (A.8) 

Evidently this system of equations has solutions for arbitrary Im(6), c, d and there- 
fore the minimal Lagrangian can be always obtained by an equivalence transfor- 
mation (|A.3| ) of fermion fields. At the quantum level this local transformation does 
not yield a nontrivial Jacobian when one applies the dimensional regularization to 
calculate it. 



Appendix B. 

In this and the following appendixes we present the results of our perturbative 
calculation. When necessary, both finite and divergent parts are given. Diagrams 
are labeled according to the figure number. 
Diagram l.a: 

-^9(A- B)A{0)U- 1 dlUU~ 1 . (B.l) 
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Fig. 2. One-loop diagrams for the vertex with one x-ficld. 



Diagram l.b: 
1 



0(A - B) [(3 + z 2 )A(0) + (1 - z 2 )A(A, B)] U^d^UU^d^UU- 1 . (B.2) 



Appendix C. 

In this appendix we show the result of the one-loop calculation of the vertex with two 
fermions and one boson line. The divergences that appear are not fully eliminated 
by a redefinition of U and additional counterterms with higher derivatives are called 
for. 

Diagram 2. a: 

- \e{A - B)A(0)U- 1 d^d 2 U)U- 1 [x^A) + x„{B) + z (x„(A) - x^B))} . (C.l) 
Diagram 2.b: 

U(A B)U- 1 { [A(0) (2(1 + z)x,{A) + (1 - z)%(B)) 

+A{A, B)(l - z 2 )x^B)] dvUU-^dpdvU 
+ [A(0) ((1 + z)%(A) + 2(1 - z)x„{B)) 

+A{A, B)(l - z 2 )x^A)] d^UU-^U} U-\ (C.2) 
Diagram 2.c: 

- 1 -6{A B) { A(0) [(1 + z) 2 (3 - z)x,(A) + (1 - z) 2 (3 + z)x,(B)] 
+A(A, B)(l - z 2 ) [x^(A)(l + z)+ - z)\ 

+ (1 - z 2 ) £ dr^(r) [A(r, B)(l - z) - A(A, r)(l + 

xU^dvUU^d^UU^dvUU- 1 . (C.3) 
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Appendix D. 

This appendix has to do with the fermion determinant and the vanishing scale 
anomaly discussed in the main body of the paper. 

There are two mutually conjugated operators in the bilinear form of the La- 
grangian (|): 



V 



(1 - z)U(r)id T + (1 + z)id T \U{t) 

(1 + Z*)U+(T)ld T + (1 - Z*)ld T (u+(t) 



(D.l) 



Therefore the total fermion determinant (the result of integration over fermions) or 
fermion loop contribution can be represented by 



(id T - ~(1 - z)UU + ){id T - - z*)UU+) 



(D.2) 



where we have restricted ourselves with unitary fields U . Now one can factorize 
out the infinite constant for free operators and find the nontrivial part in terms of 
fermion propagators: 



Z f 



|(i9 T ) 2 ||exp Tr log 1 



+ log 1- 



2id T 



(1 - z)UU~ 



2id T 

1 



{l-z*)UU A 



Cexp |-^(0)| J dr(l -z + l- z*)tr{UU+) 
Cexp(-i / dTtv{UU + ) 



(D.3) 



Herein the triangle property of the free fermion propagator has been exploited, 



(D.4) 



id T + ie 

which follows from the advanced Green function prescription by adding +ie. As 
a result, only the first order in the expansion of the logarithms in (D.3) survives 
when the functional trace operation is performed. The value 6(0) = 1/2 and the CP 
invariant choice z = —z* are employed. More rigorously this result can be obtained 
by taking the finite proper-time interval with anti-periodic boundary conditions 
for fermion fields. Then the determinant will be given by the product of discrete 
eigenvalues of the operators ( D.l ). When proceeding to the infinite line limit and 
for the advanced prescription with +ie one recovers exactly the functional presented 
in the last line of (D.3). 

Evidently, for SU(N) groups and for U(l) groups with periodic boundary con- 
ditions the exponent in (D.3) vanishes and therefore the fermion loop contribu- 
tion from the minimal Lagrangian is absent. It automatically eliminates the scale 
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Fig. 3. One-loop counterterms to the fermion propagator coming from the additional vertices. 

anomaly ~ (Lf) vac and thereby the conformal symmetry remains intact by vacuum 
polarization effects. 

For the extended Lagrangian including the higher dimensional vertices (|2^) the 
derivation of the fermion determinant is similar. The corresponding differential 
operators look as follows: 



V 



id T - -(1 - z)UU+ + -(1 - z 2 ) ( 9l - zg 2 )d v UU + d„UU + { 9l + zg 2 ) 
Z o 



-d u UU + d u UU + + 2zg z d v UU + UU + d v UU A 



U, 



& = U+ Udr - -(1 - z*)UU+ - -(1 - (z*) 2 ) [( gi - z*g 2 )Ud v U+Ud v U^ 
-(gi + z*g 2 )Ud v U + Ud v U + + 2z* g 3 Ud u U + UU + Ud u U^ 



(D.5) 

Again only the first order in the expansion of the logarithmic trace is not vanishing: 
1 



log Z f = C 



8 



(< 



drtr f 1 - z 



zg 2 d v Ud v U + + z gs ,d v UU + UU + d v UW 



+ (1 - (z*) 2 ) -z*g 2 d v Ud v U+ + z*g 3 d l ,UU + UU + d l/ UU^ 



(D.6) 



The vertices proportional g\ do not appear after tracing. Let us take the CP sym- 
metric constants z = —z*. The vertices with the coupling constant g 3 happen to be 
proportional to z + z* and therefore vanish. As to the g 2 terms they form a total 
time derivative in the CP invariant case, 

1 r°° 

logZ f = C' + - d T {l-z 2 )zg 2 d T tY{d„Ud v U + )=Q, (D.7) 

° J-oo 

the latter taking place for periodic boundary conditions. 

Thus, in spite of the fact that the higher dimensional vertices bring a scale 
dependence into the Lagrangian (the constants gj are dimensional), they do not 
generate a scale anomaly due to fermion loops iff the CP symmetry is imposed on 
the Lagrangian. 



Appendix E. 

Here we list the contribution from the additional counterterms 
propagator. 
Diagrams 3: 



|) to the fermion 



February 1, 2008 14:1 WSPC/Guidelines-IJMPA hadstr 



28 J. Alfaro, A. A. Andrianov, L. Balart, D. Espriu 



Fig. 4. One-loop diagrams for the vertex with two z-fields. 



- (A(0) - A(A, B)) [(1 + z^d^UU- 1 ^- 1 
-(l-zp-'^U-'d^UU- 1 ], (E.l) 

where 

^ = A(o)i^- (^-^d^uu^dvU - l -^d v uu~ 1 d yi d v u 

+zd„UU- 1 d^UU- 1 d u U) . (E.2) 
When retaining only the terms ~ A 2 (0) one reproduces exactly 2du . 



Appendix F. 

Here we list the diagrams that are relevant to the calculation of the one-loop vertex 
involving two fermions and two boson lines. 
Diagram 4. a: 



--6(A-B)A(0) [x„(A)x v (A)(l + + x^B)x„(B)(l - z)] U-^d^UU^.iF.l) 
Diagram 4.b: 

Its divergent and finite parts are 

±0(A - B) { A(0) [x ll {A)x v {A)2{l + z) + x^{B)x v {B){\ - zf] 
+A{A, B)x^{B)x v {B){\ - z 2 )} U^dpUU^O^dpUU- 1 . (F.2) 
Diagram 4.c: 
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Its divergent and finite parts are 

U(A - B) {A(0) [i„(A)x„(A)(l + z) 2 + x^(B)x v (B)2(l - z)\ 
+A(A,B)x tl (A)x v (A)(l - z 2 )} U^d^dudpUU^dpUU- 1 . (F.3) 
Diagram 4.d: 

Its divergent and finite parts are 

~0(A - B) {A(0) [ Xfl (A)x„(A)(l + z) 2 (3 - z) + x^(B)x„(B)(l - z) 2 (3 + z)] 
+A(A, B){1 - z 2 ) [x^A)x v (A)(l + z)+ Xfl (B)x v (B)(l - z)] 

+ (1 - z 2 ) drix^ir) + x^x„{t)) [A(t, B)(l - z) - A(A, r)(l + z)] j 

xU^dpUU^d^UU^dpUU- 1 . (FA) 
Diagram 4.e: 

Its divergent and finite parts are 

~0(A - B) {A(0) [xp{A)x v {A){\ + z) 2 (3 - z) + x„(B)x„(B)(l - z) 2 (3 + z) 

+ {1-Z) 2 (1 + Z) dT Xfl (T)x u {T) 



+A(A, B)(l - z 2 ) [xp{A)x u {B){\ + z)+ Xft (B)x v (B)(l - z)] 

+ (1 - z 2 ) £ dr [ X p( T )x„(B)A(T, B)(l -z)- ^(r)^(r)A(A, r)(l + 

xU^dpUU^d^UU^dvdpUU- 1 . (F.5) 
Diagram 4.f: 

Its divergent and finite parts are 

-\0{A - B) (A(0) [xp{A)x v {A){\ + z) 2 (3 - z) + - z) 2 (3 + z) 



-(l-z)(l + z) 2 / dr^(r)^(r) 

JB 

+A(A, B)(l - z 2 ) [xp{A)x v (A)(l + z) + x»{A)x v {B){\ - z)] 

+ (1 - z 2 ) £ dr [x^{t)x v {t)A{t, B)(l - z) - x^A)x v {t)A{A, r)(l + z)] j 

xU^dpdpUU^dvUU^dpUU- 1 . (F.6) 
Diagram 4.g: 

Its divergent and finite parts are 
\e(A - B) {A(0) [xp{A)x v (A)(l + z)(3 + z)) + x^{B)x v {B){l - z)(3 - z) 
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[A 

+ (1 - Z 2 ) I dT(x^{T)x v (T) ~ X^(t)x v (t)) 

Jb 

+2A(A,B)(1 - z 2 )x^A)x u (B)} 
xU^dpd^UU^dpdvUU- 1 . 

Diagram 4.h: 

Its divergent and finite parts are 



(F.7) 



1 



-9(A - B) (A(0) [x lx {A)x v {A)2{\ + z) 2 (3 - z) + x lx {B)x v {B)2{\ - z) 2 (3 + z)} 
+A(A, B)(l - z 1 ) [x ll {A)x v {A)2{\ + z) + i„(B)i„(B)(l - zf 

+Xp{A)x 1/ {B){l- z 2 )-{\- z 2 ) / dTX^{ T )x u {T) 

JB 

f A 

+ (1 - z 2 ) / dr [A(r, B) (i„(T)£„(T)2(l - *0 

JB 

+^(r)x,(r)(l - z 2 ) + i„(T)i v (B)(l - z) 2 ) 

-A(A, t) (i M (T)i„(T)(l - z 2 ) + ^(r)x,(r)2(l + z) + i M (A)i„(T)(l + z) 2 )] 

-(l-Z 2 ) 2 J dn J dT2X ll (T 1 )x 1/ (T 2 )A(T 1 ,T 2 )^ 

xU^dpUU^d^UU^dMU^dpUU- 1 . (F.8) 

All divergences are removed by combining the rcnormalization of U and the 
additional counterterms determined from the one loop vertex with one external 
boson lines. No new counterterms are required. 



Appendix G. 

Two loop diagrams for the fermion propagator. Only the divergent parts are needed. 
Diagram 5. a: 

Its contribution is fully divergent, 

- h(A - B)A 2 (0)U- 1 dld 2 p UU- 1 . (G.l) 

Diagram 5.b: 

Its divergent parts are 

^6{A - B) {A 2 (0)(3 + z 2 ) + A(A, B)A(0)(1 - z 2 )} U^dldpUU^dpUU^ 1 . (G.2) 

Diagram 5.c: 

Its divergent parts are 

^9(A - B) {A 2 (0)(3 + z 2 ) + A(A, B)A(0)(1 - z 2 )} U^dpUU^dldpUU' 1 . (G.3) 
Diagram 5.d: 
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(b) {c> 



W (f) 



00 0) 



Fig. 5. Two-loop diagrams for the propagator. 

Its divergent parts are 

~0(A - B) { A 2 (0)(3 + z 2 ) + A(A, B)A(0)(1 - z 2 )} 
xU^dpUU^dlUU^dpUU- 1 . (G.4) 
Diagram 5.e: 

Its divergent and finite parts are 

-^6»(A - B) {A 2 (0)(11 + z + hz 2 - z 3 ) + A(A, B)A(0)2(1 - z 2 )(l - z) 
+A 2 (A,B)(l- z 2 )(3 + z)} U^dpUU^dpUU^d^dpUU- 1 . (G.5) 
Diagram 5.f: 

Its divergent and finite parts are 

~0(A - B) {A 2 (0)(11 - z + hz 2 + z 3 ) + A(A, B)A(0)2(1 - z 2 )(l + z) 
+A 2 {A, B){1 - z 2 )(3 - z)} U^dpdpUU^d^UU^dpUU- 1 . (G.6) 
Diagram 5.g: 

Its divergent and finite parts are 
\e{A - B) {A 2 (0)(3 + z 2 ) + A 2 (A, B)(l - z 2 )} U^dpdpUU^dpdpUU- 1 . (G.7) 

5 

Diagram 5.h: 
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Its divergent and finite parts are 
±8(A-B){A 2 (0)(3 + z 2 ] 
+(l-z 2 )(5-z 2 )A 2 (i,B) 



^-9(A - B) { A 2 (0)(3 + z 2 )(7 + z 2 ) + 2(1 - z 2 )(3 + z 2 )A(A, B)A(O) 



+2(1 - z 2 ) 2 jf dn J J dT 2 d Tl A{n,T 2 )d T2 A(T 1 ,T 2 )^ 

xU^dpUU^d^UU^d^UU^dpUU- 1 . (G.8) 
Diagram 5.i 

Its divergent and finite parts are 

~9{A - B) [A 2 (0)(1 + z 2 ) + A(0)A(A, B)(l - z 2 )] 
xU^dpUU^dpdpUU^dpUU" 1 . (G.9) 
Diagram 5.j: 

Its divergent and finite parts are 

-}-9(A - B) { A 2 (0)(1 + z 2 )(9 - z 2 ) + 2(1 - z 2 )(3 - z 2 )A(A, B)A{0) 
lb 

+ (1 - z 4 )A 2 (A, B) - (1 - z 2 ) 2 / dr (A(A, r) - A(A B)) A(r, B) 

Jb 



J dn J dr 2 9 T1 A(Ti,T2)a r2 A(ri,T 2 )| 



xU^dpUU^dpUU^dpUU^dpUlf- 1 . (G.10) 

The first integral is assembled into a non-singular expression. The second integral 
is prepared to have a weaker singularity than A 2 (0) and it reveals a singularity of 
A(0)-type. 

Appendix H. 

Let us first understand the singularities in J (A, B) with the help of Euclidean 2-dim 
scalar field propagator (051) in the cutoff regularization, 



A(t.o-) = -a'ln[(r 2 + a 2 + R 2 )fi 2 ] , R^O. (H.l) 

The regularization smears the log singularity, the normalization is taken to provide 
-<9 2 A(t, a) = iwa'S(r)S(a) for R -> 0. Thus on the boundary, 

A(r,CT = 0) = -a'ln[(r 2 +i? 2 ) i u 2 ] , d a A(r, a = 0) = 0, (H.2) 

the latter being in accordance with Neumann boundary conditions for open strings. 
Evidently, the relation between divergences in DR and the present regularization 
is: 

A(0) ~ -a' In [iiVl *->- — . (H.3) 
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Of course this propagator is not necessarily exact in its finite part. But the diver- 
gences extracted with its help should be universal. 
With this ansatz it is easy to find that 

J {A — B) = a' 2 j-"^— - - 2 In (i?V) j + regular terms. (H.4) 

The first divergence is power-like and we neglect it (it should not appear in DR 
calculation). The second term is logarithmically divergent ~ 2a'A(0) and we retain 
only this one. 

Now let us perform the same calculation in the Dimensional Regularization. 
First we define the integral ([39]) in 2 + e dimensions. Evidently we can do it con- 
sistently for the second expression in (|39|). Namely, we replace t — * |r| = t with f 
being a 1 + e dimensional vector and integrate over the sphere |f | < (A — B), 



1 



2 



J e (A-B) = -- d^r^(A-B-\f\) A(|f|) . (H.5) 

J\t\<(A-B) l j 

Next we insert in (H.5) the derivative of the string propagator jl5| ) in 2 + e dimen- 
sions, 

A'(*) = - a 'er(!U^VV- (H.6) 



2J \ ip 



which leads to 



2 V ' \ \2)J ^ V V 



-2( 



rA-B 

/ dt(tfi) £ {A - B -t)t- 2 - t , (H.7) 
Jo 



where the angular volume, 



Vi+t - (H- 



After integration one finds the following expression 

2 ((A - B)fiy/n)~ 



UA-B) = -(aV £ ) 2 ( e r(f)) 



,o 4(aV~ £ ) 2 



Thus we have reproduced the same value of the constant cy. 



m e(l + e) 
2aV" e A(0). (H.9) 



Appendix I. 

Let us take the most general set of operators which can appear in the equations of 
motion (E.o.M.) with arbitrary constants. The equations of motion of dimension 
two ( |2(i| ) are assumed to hold and therefore we do not include any vertices containing 
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the D'Alambertian d 2 . We would like to find out a set of constants which supports 
the unitarity relation (|2l|), USW = —SlfW. The results are presented in the 
Table 4. 



Table 4. Comparison between the coefficients of the different chiral field 
structures in the equations of motion, their Hermitean conjugates and 
those ones derived from a local Lagrangian (142) 





CF structure 




E.o.M. 


(E.o.M.) 1 " 


X-lagr. 




p p up 




01 


-02 


-2(2^1+^2) 




pp p p 




ai 


-ai 


-2{2K X + K 2 ) 




p pp p 




03 


-a 3 


-4K 2 


P 


P P 


/' 


«4 


ai + a 2 + a4 


2[(1 - 2 2 )^i + K 2 ] 


P 


P P 


P 


OS 


<23 + «5 


-2z 2 K 2 


P 


P P 


P 


«(i 


Cll + Ou-2 + «3 + a 6 


4{Ki + K 2 ] 



One can see that the unitarity of 5U^ is provided for z 2 = — 1 only if 
a\ = a 2 = -(4Ki + 2K 2 ), a 3 = -4K 2 , 04 = - -jiflx + a 2 ), 

a 5 = -ia 3 , a 6 = -i(ai + a 2 + a 3 ). (LI) 

Thus unitarity is achieved when the operators in the equation of motion are derived 
from the local Lagrangian ( f4l| ) and vice versa. 

Appendix J. 

In this appendix we shall explore a somewhat different approach and see how the 
results are fully equivalent to those presented in the text. We shall integrate out 
the 'quarks' and derive an effective action in terms of external sources. 

To derive the effective action let us supplement the Lagrangian (g) with the 
external sources for fermion fields 

Lf = Lf + JliPr + Jripl + ipL Jr + i>R Jl- (J-l) 
As this Lagrangian is quadratic in fields the effective action, 



is supported by the solutions of classical equations, 



ip R +~(l + z)U- 1 Uip R =iU- 1 J R , (J.3) 

i>R + \{l + z^rII^U^ 1 = -iJ R tf~\ (J.4) 

i>L + 5(1 - z^U^uHl = iU ] ~\j L , (J.5) 

Ql + ^-z^lUU" 1 = -U L U- 1 . (J.6) 
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The solutions read 



f dn 


Texp (-£ 


J — oo 




xiU-HrJMn), 


/ dn 

J — oo 


Tmr {-L 



dT2^(l + z)U- 1 (T 2 )U(T 2 )J 



dr 2 i(l-z*)£/t 1 (r 2 )C/t( T2 ) 



(J.7) 



(J.8) 



and their complex conjugated partners. 

Therefore the effective action takes the simple form: 



dTJ L {r) / dri 



Trxp ( - / dT2~(l + z)U- 1 {T 2 )U(T 2 ) 



In turn the full fermion propagator takes the form 



(J.9) 



1 



1 



(1-z) drUU-^x^T)] 



9(ti - t 2 ) 



~{l + z) / drU-'Ulx^T) 



Texp 



(J.10) 



We stress that two pieces with T-exponentials in the last equality are identical. 
Let us expand the first expression for the propagator in ( J.lOj ) 



1 + \{1-Z) £ dTUU-^X^T)] 



1 



+ -(1-Z) 2 / dT dT'UU- 1 [x^T)]UU- 1 [x^T')] + 



(J.ll) 



The second order in the expansion is sufficient to analyze one-loop divergences 
of the propagator. In turn we can develop the perturbative expansion around a 
background, x(t) = xq + X ( T ), 



+^x li (T 1 )x u (n)[-(d tt d u U)U- 1 (x ) + 2{d ti U)U-\d„U)U- L (x )j. (J.12) 

Evidently the term with two derivatives on U comes out only from this part of 
expansion. Further on one should also evaluate: 



UU-^x^t)} = x^UU- 1 {x )+x^ [d^UU-Xxo) 
-d„UU- l d v UU- x {x )] +••• 



(J.13) 
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Now we insert the above two expansions into Eq. (J.ll) and retain only terms 
quadratic in 

U-^xo^i-d^UU-^xo) + 2d^UU- 1 d v UU- 1 (x )]Z tl {T 1 )x v { Tl ) 
z)U- 1 d tl UU- 1 d v UU- 1 {x Q ) f 1 dTx^n)xv(T) 

- z)U- 1 [d ll d v UU-\x ) - d^UU^dvUU- 1 ^)] I 1 (It^(t)x v (t) 



hi - zfU^d^UU^dvUU-^xo) [ 1 dr [ dT'^(r)i v (T') 



(J.14) 



After integration over r, r and averaging in i M (t) with the help of formulas for the 
string propagator, having in mind that A(0) = (i.e. the contribution of the third 
line is equal zero) one obtains the 1-loop part in the form written in the main text, 

1 



(ri - t 2 )U- x 



A(0) \-diU 



+ ^^d ti UU- 1 d^u\ + (1 Z2) A(t! - T 2 )d li UU- l d»U 



If- 1 . (J.15) 



